Subgap resonances and conduction channels in mesoscopic superconducting devices. 
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Peaks associated to multiple Andreev reflections have been used to identify the number of conduction 
channels, and their transmitivity, in mesoscopic superconducting junctions. We analyze the influence 
of the detailed shape of each channel, and the number of weak links within it, on the final results. 
Connexions with the statistical interpretation of conductance histograms is also made. 



PACS numbers: 61.16-Gh, 62.20.Fe, 73.40.Cg 

Many experiments show that the conductance of nar- 
row constrictions between different metals is quantized 
jfl-jlj). The nature of the structures which give rise to 
this effect is not well known. The forces transmitted 
through the constriction show discontinous jumps cor- 
related with the conductance steps fH| , suggesting that 
abrupt atomic rearrangements give rise to the conduc- 
tance steps. Direct observations show irregular, and even 
amorphous structures near the contact |13| [L{M . 

The simplest explanation of the observed quantization 
describes the region near the constriction in terms of a 
fixed number of almost fully open channels. Such a pic- 
ture is consistent with models in which the constriction 
is described in terms of nearly free electrons |2(i|-|25[| . Re- 
lated schemes have been used to describe the force oscil- 
lations observed experimentally |2^-|29[| . Statistical aver- 
ages based on models where the connexion between the 
channels within the constriction and the electrodes are 
allowed to vary randomly, support this picture A 
model of an n-channel constriction coupled to TV-channel 
electrodes, shows the desired features, when n <C N, and 
the reflection and transmission of electrons from the elec- 
trodes are taken to be random matrices fl3l| . 

An alternative picture, based on local properties of the 
atoms at the constriction has also been put forward [§2| . 
The conductance is defined using localized wavefunctions 
at the atoms within the constraint. When the constric- 
tion narrows to a single atom at some point, conduc- 
tances near a single quantum are obtained. The bands 
of metals such as Pb or Al arise, basically, from the hy- 
bridization of s and p atomic orbitals, while those of Na 
or Au are mainly derived from s orbitals. Hence, the 
number of channels through a one atom constriction in 
a nearest neighbor tight binding scheme is at least three 
for Al or Pb, while it is one for Au. 

These two pictures can be adapted to give similar re- 
sults for the total conductance of the constriction, at 
least when a single quantum is observed. The properties 
of superconducting constrictions, however, are expected 
to differ. The FV characteristics of a superconducting 
constriction described by a single channel with one bar- 
rier are significantly different from those of three or more 
channels, even if the total conductance in the normal 
state is the same. The Andreev spectrum below the su- 
perconducting gap |33] has a much richer structure in the 



second case [[Mj. Experimental results show, very conclu- 
sively, that the Andreev spectrum in Al and Pb is better 
fitted by a model of three or more superconducting chan- 
nels in parallel, with a suitably adjusted transmission co- 
efficient |m| . Current investigations try to elucidate if a 
similar distinction can be made in other quantities which 
do not depend linearly on the total conductance, such as 
the critical current or the noise spectrum. 

We will now focus on the difference between the two 
models described above for the case when a single con- 
ductance quantum is observed experimentally. A sta- 
tistical description based on a nearly free electron pic- 
ture assumes that the constriction contains a single elec- 
tronic channel. Because of the properties of the (ran- 
dom) connexions to the bulk electrodes, in most cases, 
the transmission of this channel is near one. The larger 
the number of channels in the external electrodes, the 
more likely it is to find a combination of the channels 
in the electrodes with a good matching to the channel 
in the constriction, explaining the observed clustering of 
conductances near Making a transformation to a 

basis in which these combinations are singled out, the 
model can be reduced to a one dimensional single chan- 
nel problem, without loss of generality. The model, how- 
ever, requires two barriers, at the two connexions of the 
weak link with the electrodes, as shown in fig.(|l|). On 
the other hand, a model based on atomic orbitals with 
nearest neighbor hybridization, in which the constraint 
narrows to a single atom, can be made equivalent, in the 
normal state, to a number of electronic channels of the 
order of the number of orbitals in one atom. Because of 
the small size of the constriction, the description of the 
transmission through each of these channels by a single 
barrier is highly plausible. Note, however, that the for- 
mation of the superconducting condensate can introduce 
changes in this picture at scales comparable to the gap. 

In the following, we analyze the I-V characteristics of 
superconducting structures made up of a single channel 
with two barriers, which separate the constriction from 
the bulk electrodes. We can always hybridize the chan- 
nels in the bulk electrodes (regions I and III in figure ([j])) 
in such a way that only one is coupled to the channel 
within the constriction. Thus, without loss of general- 
ity, the problem is reduced to that of a single channel 
interrupted by two barriers. 
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FIG. 1. Sketch of the representation of the constraint ana- 
lyzed in the text. Channels in regions I, II and III are assumed 
to have perfect transmission. The complexity of the junction 
is represented by the random matrices connecting the differ- 
ent regions. 

We use the scattering formalism presented in p5[ , 
which is equivalent to the Green's function method for- 
mulated in p6| and used in the interpretation of experi- 
mental I-V curves . Each barrier is defined by a scattering 
matrix: 



Si = 



t* 



(1) 



where r and t are complex numbers which satisfy \r\ 2 + 
\t\ 2 = 1, and which are denned up to a global phase. 
In ([!]), we assume time reversal symmetry. The elec- 
tronic wavefunctions at the channels in the electrodes and 
within the constriction, fig. ([!]), are described as plane 
waves, ^k(x) ~ e tkx . The quasiparticle energies, in 
the normal state, are linearized around the Fermi level, 
£fc = £f + vp(k — hp), where cf is the Fermi energy, vf 
is the Fermi velocity, and kp is the Fermi wavevector. 
Each wavefunction has different phases at the two barri- 
ers. The overall phaseshift includes a contribution from 
the region described in terms of perfect channels, and 
a contribution from the matching zones, described phc- 
nomenologically by the scattering matrices. The phase- 
shift can be written as an energy independent part plus 
an energy dependent contribution, which we parametrize 
as: 



(e) = /- 



(2) 



where / is dimensionless, and A is the superconducting 
gap. If the transmission through the constraint is bal- 
listic, we have / = where L is the length of the 
region between the barriers. A more general situation is 
discussed later. 



The model has four adjustable parameters: the abso- 
lute values, \ti\ 2 , the phase between n and which we 
label 4>o, and / as defined above. 

In the normal state, the system formed by the two 
barriers in series has an energy dependent transmission: 



T(e) = 



T t T 2 



1 + RxR 2 - 2/R^cos (0o + 2/4 



(3) 



where Ri = \ri\ 2 and T$ — \ti\ 2 are the reflection and 
transmission coefficients at each barrier. 

The most general form of the wavefunctions at regions 
I, II and III, is: 
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where V\ and Vi are the voltage drops at the two bar- 
riers, a™(e) = a(e + 2meVi + ZneV^) is the Andreev 
reflection amplitude, and Jo i n the first equation gives 
the incoming current. The matching conditions lead 
to a set of matrix equations between the coefficients 
A, B,C, D, E, E, F and F in the wavefunctions. It can 
be shown, after some algebra, that the only finite coeffi- 
cients have m = n. This implies that physical quantities 
depend on the total voltage drop, V\ + V 2 - This is a 
consequence of our neglect of inelastic scattering in re- 
gion II. Finally, the coefficients in region II, E, F, E and 
F, can be written in terms of the other four, leaving a 
set of equations formally equivalent to those of a single 
junction. They can be solved using recurrent fractions, 
as described in |5?J . 

The main difference with the case of a single barrier is 
the energy dependence of the phase, parametrized by /. 
Multiple Andreev scattering is determined by the trans- 
mission coefficient of the constriction at different ener- 
gies. Hence, it is reasonable to expect that this energy 
dependence has a significant influence on the I-V charac- 
teristics below the gap energy. 
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FIG. 2. I-V characteristics at low voltages of supercon- 
ducting junctions. Full curves are obtained with the parame- 
ters used in[34]. Upper curve: two channels of transmissions 
T\ = 0.74 and T2 = 0.11. Lower curve: three channels of 
transmissions Ti = 0.46, T 2 = 0.35 and T 3 = 0.07. Squares 
and circles are approximations obtained with the double bar- 
rier model described in the text. The broken curves are results 
for a single channel, one barrier model with the same normal 
state conductance. 
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FIG. 3. Josephson current, as function of the phase differ- 
ence between the superconducting electrodes, of the two sets 
of barriers described in fig. (2). Full line: single channel with 
one barrier and the same transmission coefficient as the set of 
parallel channels considered in[34]. Broken line: Results for 
two or three channels in parallel. Dash-dotted line: Results 
for a single channel with two barriers in series. The upper 
panel corresponds to the two channel case of fig. (2). The 
lower panel corresponds to the three channel case. 



In fig.(g), we show the best fittings that we obtain to 
the I-V characteristics analyzed in ( J|4|]) in terms of dif- 
ferent channels in parallel. The method of minimization 
of the mean square deviations that we use leads to more 
than one possible fitting of similar quality. The param- 
eter / is 0.08 for the upper curve (which corresponds to 
a two channel model in p|) and 0.1 for the lower curve 
(which corresponds to a three channel model in 



Another quantity which has a non linear dependence 
on the transmission, and which can be used to elucidate 
the detailed structure of the constraint, is the critical cur- 
rent. In fig-© we show the dependence of the Josephson 
current [|38f| on the superconducting phase between the 
electrodes for the two cases analyzed in fig.(^) and for a 
single channel and one barrier of the same normal state 
conductance. The maximum value gives the critical cur- 
rent of the junction. Our results deviate from the single 
channel one barrier case in the opposite direction to those 
obtained by combining different channels. 



In order to relate / to the structure of the constraint, 
we assume that the motion of the electrons in the con- 
tact zones is diffusive, with a mean free path, The 
dependence of the phaseshift on distance and frequency 

T 2 

goes as where D 
Replacing w by A, we obtain: 



^ is the diffussion coefficient. 



AL 

TivfI 



(5) 



where L is the length of the constriction, and £o is the 
coherence length at zero temperature. Taking £o ~ 10 4 A 
for Al, the values of / that we obtain, ~ 0.1, suggest 
values of the order of L ~ 300 A and I ~ 100 A. These 
numbers are consistent with experimental observations 

The enhancement of the Andreev reflections that we 
find in a double barrier geometry is similar to the results 
obtained for transport through a single level between su- 
perconducting electrodes Q . 

In summary, we find that existing experiments cannot 
rule out a description of many of the weak links which 
show conductance quantization in terms of a small num- 
ber of channels (one for each quantum) randomly coupled 
to larger systems described by a much higher number of 
channels. This description is consistent with the observed 
statistical distribution of the normal state conductance. 
Our fitting of experimental results require constrictions 
of typical dimensions L ~ 300 A, and strong disorder, 
with mean free paths I ~ 100 A. The size of the constric- 
tion, and the existence of strong disorder, are compatible 
with existing experiments. 
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The model analyzed here, however, needs not apply 
to all observed cases, as it only attempts to describe an 
average situation. When the constriction size is of the or- 
der of a single atom, it is plausible conducting evanescent 
waves (in our language) play some role. If their contri- 
bution to the conductance is a significant fraction of the 
total value, a description based on extended wavefunc- 
tions is almost indistinguishable from one based on lo- 
calized atomic orbitals. On the other hand, the existence 
of well defined conductance steps near integer multiples 

2 

of %- seems to us more consistent with a dominance of 
situations with almost fully conducting channels within 
the constriction. 
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